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TOMOSYNTHESIS IMAGING SYSTEM AND METHOD 

RELATED APPLICATIONS 

This application claims priority to provisional application serial number 60/446,784, filed 
5 February 12, 2003 and entitled Tomosynthesis Imaging System and Method, which is 
incorporated by reference herein. 

BACKGROUND OF THE INVENTION 
1 . Technical Field of the Invention 
10 The present invention relates to a system and method for imaging a target element using 

tomosynthesis. More specifically, the invention relates to a system, method and computer 
program product for creating a three-dimensional image of target elements from a plurality of 
radiation absorbance projection images taken from different angles. 

15 2. Background 

Imaging of a patient's tissue has become a common screening and/or diagnostic tool in 
modern medicine. One example of such imaging is mammography, or the imaging of a patient's 
breast tissue. Breast cancer remains the most common cancer among women today, however, at 
this time there is no certain way to prevent breast cancer and the best strategy for dealing with 

20 breast cancer is early detection of the cancer so that it may be treated prior to metastatic spread. 
Accordingly, it is important for patients to have access to imaging techniques and systems that 
will detect very small cancers as early in their development as possible. 

Conventional mammography involves an x-ray examination of the breast, typically using 
25 a fluorescent panel that converts transmission x-rays from a breast into visible light photons that 
expose a film. While screening using conventional mammography has been shown to reduce 
breast cancer deaths by approximately 30 to 50%, this imaging technique lacks the dynamic 
range that would allow it to detect small or hidden cancers, and thus permit therapy that can 
improve survival rates further. In particular, conventional mammography techniques suffer from 
30 the limitation that three-dimensional anatomical information is projected onto a two dimensional 
image. Because of this, "structure noise" such as overlapping breast tissues makes it difficult to 
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perceive and characterize small lesions. This can result in a 10 to 30 % false-negative diagnosis 
rate, especially where the cancer is masked by overlying dense fibroglandular tissue. 

A three-dimensional approach to imaging could allow for the separation of overlying 
5 tissue and thus improve correct diagnosis rates for diseases such as breast cancer; however, 
three-dimensional imaging has not yet been applied for this purpose in the general population. 
The most widely used three-dimensional x-ray imaging technique is computed tomography 
("CT"). A CT scanner contains a rotating frame that has one or more x-ray tubes mounted on 
one side and one or more detectors on the opposite side. As the rotating frame spins both the x- 

10 ray tube and the detector around the patient, numerous projections of the x-ray beam attenuated 
by a cross section slice of the body are acquired. These projections are then used to reconstruct 
cross-sectional images of the body. Despite the fact that CT has been found useful in detecting 
lesions in the breast, it is not suitable as a technique for regular breast imaging due to the high 
dose required to take a number of projections (approximately 100 to 1,000 projections) and the 

1 5 low spatial resolution (on the order of a millimeter). In addition, the CT projections mix 

attenuation effects from other organs of the body (such as those within the chest cavity) with the 
attenuation of the breast, which can distort information about the breast and causes these 
interposed organs to be irradiated. Still further, the cost of CT scanning is too high to permit its 
use as part of an annual exam. 

20 

A three-dimensional imaging approach called "tomosynthesis" has also been developed. 
Tomosynthesis is a technique that allows the reconstruction of tomographic planes on the basis 
of the information contained in a series of projections acquired from a series of viewpoints about 
the target object. They need not be regularly spaced, numerous, or arranged in any regular 
25 geometry. The tomosythesis technique is promising in that it may provide improved spatial 
differentiation of nearby tissues at very high resolution comparable to projection 2D imaging, 
with limited radiation. The problem of 3D reconstruction from tomosynthesis projections has 
been described as intractable by those skilled in the art. 

30 In order for a three-dimensional imaging technique to be successful in medical diagnosis 

and other applications, it should offer: 
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• Sufficient spatial resolution and contrast resolution to detect and characterize, for 
example, breast cancers; 

• Minimum radiation dose to a patient; 

• Fast image acquisition; 
5 • Cost effectiveness; and 

• 3D reconstruction that can be performed effectively. 
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SUMMARY OF THE INVENTION 

In one aspect, the invention provides a method that enables to the use of tomosynthesis to 
efficiently provide accurate three-dimensional imaging of a target element. This method 
involves acquiring radiation absorbance images of the target element through a limited plurality 
5 of angles and applying an iterative reconstruction algorithm to generate the three-dimensional 
reconstruction of the target element. This method can gain further accuracy where the iterative 
reconstruction algorithm is applied using cone beam forward projection and back projection. 

In a further aspect of the invention, a system for three-dimensional tomosynthesis 
10 imaging of a target element is provided having an image acquisition element and a processor. 

The image acquisition element obtains a plurality of images of the target element from a plurality 
of angles and includes a radiation source that is positionable at a plurality of angles with respect 
to the target element and a radiation detector. The radiation detector is positioned so as to detect 
radiation emitted by the radiation source passing through the target element and determine a 
15 plurality of attenuation values for radiation passing through the target element to establish a 
radiation absorbance projection image of the target element for a particular radiation source 
angle. The processor is configured to apply an iterative reconstruction algorithm to the radiation 
absorbance projection images of the target element obtained from a plurality of radiation source 
angles to generate a three-dimensional reconstruction of the target element. Again, the system 
20 can gain further accuracy where the iterative reconstruction algorithm is applied using cone- 
beam forward projection and back projection. 

In a still further aspect of the invention, a computer program for three-dimensional 
tomosynthesis imaging of a target element is provided. The three-dimensional images are 

25 created from a plurality of radiation absorbance projection images obtained at different angles 
from an image acquisition element having a radiation source positionable at a plurality of angles 
with respect to the target element and a radiation detector. The radiation detector is positioned 
so as to detect radiation emitted by the radiation source passing through the target element and 
determine a plurality of attenuation values for radiation passing through the target element to 

30 establish a radiation absorbance projection image of the target element for a particular radiation 
source angle. The computer program code is embodied in a computer readable medium and 
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includes computer program code for applying an iterative reconstruction algorithm to the 
radiation absorbance projection images of the target element obtained from a plurality of 
radiation source angles to generate the three-dimensional reconstruction of the target element 
wherein the iterative reconstruction algorithm is applied using cone-beam forward projection and 
5 back projection. 

In specific embodiments of any of these aspects of the invention, the radiation absorbance 
images can be acquired by transmitting x-ray energy from an x-ray source through the target 
element to an x-ray detector and the x-ray detector may have a plurality of detector pixels. The 

10 three-dimensional reconstruction of the target element may be represented as an array of voxels 
having a uniform or non-uniform size in three-dimensions. The forward projection step may 
then be implemented by ray tracing from the x-ray source to a detector pixel and the forward 
projection of the target element is obtained by repeating the ray tracing for each detector pixel to 
calculate an attenuation value for each voxel. The back projection step can be implemented by 

1 5 locating detector pixels containing attenuation information relating to a selected voxel and using 
those pixels to update the attenuation value of the selected voxel. The back projection step can 
further include performing a back projection for at least each voxel corresponding to a three- 
dimensional reconstruction of the target element. In the enumerated aspects of the invention or 
in any of their embodiments, the iterative reconstruction algorithm may be a maximum 

20 likelihood algorithm and the maximum likelihood estimation can be implemented using an 
expectation-maximization algorithm. 

The invention is particularly useful for creating three-dimensional reconstructions of 
animal and more particularly human tissue. In one preferred embodiment, the invention is 
25 employed in mammography to create a three-dimensional reconstruction of the breast tissue of a 
human female patient. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The invention will be more fully understood from the following detailed description 
taken in conjunction with the accompanying drawings: 

FIG. 1 is a diagram illustrating the geometry of a tomosynthesis system of the invention; 

5 

FIG. 2 is a top view of the coordinate system of a tomosynthesis system of the invention; 

FIG. 3 illustrates a forward projection of the tomosynthesis system of the invention; 

10 FIG. 4 illustrates the path length of an x-ray beam in a voxel in the tomosynthesis system 

of the invention; 

FIG. 5 illustrates a projection of the path length of FIG. 4; 

15 FIG. 6 illustrates an exception to the projection of FIG. 5; 

FIG. 7 illustrates a back-projection step of the invention; 

FIG. 8 A illustrates a phantom used to test a system of the invention; 

FIG. 8B illustrates a feature plate that makes up a portion of the phantom of FIG. 8A; 



20 



FIGS. 9 A, 9B, and 9C illustrate structural noise reduction in (A) a projection of an ACR 
phantom; (B) a projection of a mastectomy specimen/ACR phantom; and (C) a reconstructed 
25 ACR phantom feature layer; 

FIG. 10 is a film-screen mammogram of a patient's tissue; and 

FIGS. 1 1 A, 1 IB, and 1 1C illustrate slices of a reconstructed volume of the same tissue at 
30 three different depths. 
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DETAILED DESCRIPTION OF THE INVENTION 

The systems and methods of the present invention address the needs of the art by 
providing tomosynthesis apparatus and techniques for imaging target elements that overcome the 
problems of conventional three-dimensional imaging systems. The present invention enables 
5 the use of tomosynthesis to efficiently provide accurate three-dimensional imaging of a target 
element in which overlapping sub-elements having differing attenuation characteristics by 
applying a 3D reconstruction algorithm having a novel combination of features. The algorithm 
can employ a cone-beam geometry lacking in geometric simplification such as parallel-beam 
based approximation methods. The algorithm can further apply the cone-beam geometry in an 

10 iterative forward-projection and back-projection method based on maximum-likelihood image 
estimation using an estimation-maximization algorithm. The invention is applied below to one 
preferred embodiment in which the system is used for tomosynthesis mammography; however, 
the invention will be useful in a variety of three-dimensional imaging situations. For example, 
the invention can be applied to a variety of patient imaging problems such as heart imaging, or 

15 imaging of the soft tissues or bones of the hand. The imaging system of the invention can be 
used for diagnoses (as is described below for tomosynthesis mammography) or it may be used 
for other applications such as three-dimensional modeling for the purpose of fitting an implant 
(whether orthopedic, such as a hip or knee implant, an artificial heart, or other type of implant) or 
for use in surgical navigation systems. What follows is a description of one preferred 

20 embodiment of the invention. 

1. Tomosynthesis Mammography System 

Tomosynthesis mammography is a three-dimensional breast imaging technique. It 
involves acquiring projection images of a breast at a plurality of viewpoints, typically over an arc 

25 or linear path. Three-dimensional distribution of x-ray attenuation coefficient of the breast 

volume is reconstructed from these projections. A prototype tomosynthesis system 10 for breast 
imaging is illustrated in Figure 1. In this exemplary system, eleven projections are acquired by 
moving the x-ray tube 12 over a 50° arc (-25° to +25°) above the target element, in this case 
breast tissue 18 which may be compressed by compression paddle 16, in 5° angular steps about 

30 axis of rotation 14. Breast tissue 18 and digital detector 20 are stationary during the image 
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acquisition. Certain characteristics of this exemplary embodiment of a tomosynthesis system of 
the invention are described below: 



Spatial resolution and contrast resolution: The tomosynthesis system uses an amorphous- 
5 Silicon-based flat panel digital detector 20 on which a Csl crystal phosphor is grown epitaxially. 
It reads out 2304X1800 pixels (lOOjam pixel pitch) via a TFT array. The detector has a linear 
response over exposure levels up to 4000mR and 12 bits of working dynamic range. Each plane 
of the 3D reconstruction has about the same resolution as the detector (lOOum) but the depth 
resolution is on the order of a millimeter. 

10 

Dose: The target/filter combination is Rh/Rh and the accelerating potential is 25~33kVp to 
image breasts with 3~8cm range of thickness. The total x-ray dose for acquiring 1 1 projections 
is approximately 1.5 times of that used for one film-screen mammogram. Each projection is a 
low dose breast image (approximately 1/1 1 of the does per projection). 

15 

Patient motion: Patient motion is reduced by fast image acquisition. Using cone-beam x-ray 
geometry and area detector, a projection of the whole breast can be recorded with one x-ray 
exposure at each angle. For each projection, the exposure time is 0.1~0.2s and detector readout 
time is about 0.3s. Rotation to the next angle is performed during the detector readout. The total 
20 image acquisition time for 1 1 projections is about 7sec. Breast compression also helps to reduce 
patient motion. 

Image acquisition geometry: The design of the tomosynthesis system can be based on the 
conventional mammography system. The MLO views have been used in most cases since it 
25 provides the most complete coverage of the whole breast. 

2. 3D Reconstruction Algorithm 

Tomosynthesis can take advantage of the high efficiency of a digital detector in acquiring 
low dose breast images. Prior to the present invention, appropriate reconstruction methods that 
30 make good use of the low dose projections and the acquisition geometry of the tomosynthesis 

system 10 have not been deployed. For an initial evaluation, Niklason implemented a "shift-and- 
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add" method that is similar to backprojection [Niklason et al, 1997]. Methods used by others 
[Chakraborty et al, 1984; Haaker et al, 1985; Suryanarayanan et al, 2000] essentially did not 
handle the limited statistics in low dose projection images. In theory, they were not suitable in 
the case of limited number of projections and limited angular range. Therefore, the three- 
5 dimensional information extracted by these methods was limited, which resulted in poor quality 
reconstructions. 

The Maximum Likelihood (ML) algorithm is an iterative reconstruction method [Rockmore, 
1977; Shepp et al, 1982; Levitan et al, 1987; Herbert et al, 1989; Browne et al, 1992; Manglos et 
10 al, 1995; Pan et al 1997; Zhou et al, 1997]. It is well suited for tomosynthesis reconstruction, 
which is an ill-conditioned problem (only 1 1 low dose projections are available). The ML 
algorithm incorporates the stochastic nature of the x-ray transmission process so that the 
statistical noise in projection images is taken into consideration in the case of low x-ray flux. It 
also incorporates the information of the object into the reconstruction in the form of constraints. 



In ML reconstruction, the Likelihood function, which is the probability of obtaining the 
projections Y obtained in a measurement, given a certain model for the three-dimensional map 
of attenuation coefficients u is: 



20 The ML solution is the 3D reconstruction that maximizes the probability of the measured 

projections. Because an analytical solution is usually intractable, an iterative algorithm is a better 
choice. The incident and transmitted x-rays follow Poisson statistics and the log-likelihood is 
described by: 



25 where u is the linear attenuation coefficient; N s is the number of incident x-ray photons to 

projection pixel i, before attenuation; Y t is the number of transmitted x-ray photons to projection 
pixel i, after attenuation; l (j is the path length of beam ray i in the object (reconstruction voxel j; 



L = P(Y\u) 



(1) 



LnL = Yu(r N i e ~ <lM>i ~ Y i < ! > u >« +^ lnN i ~ lnY i) 



(2) 



and </,w > i 



, = ^lyUj is the total attenuation along beam ray to pixel i. 
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The algorithm by Lange and Fessler [Lange and Fessler, 1995] can be selected to solve 
the ML problem. At the n-th iteration, the value of an object voxel \x is updated by: 

H< n+,) = uf + Uf — — — (3) 

i 

5 where the notations are the same as above. 

3. Implementation of the Cone-beam Reconstruction 

Cone-beam forward projection and back projection can form the basis for iterative 
reconstruction according to the invention. At the forward projection step, the projection images 
10 at 1 1 angles are calculated based on the current 3D reconstruction model. At the backprojection 
step, the calculated projections and the measured projections are compared and the 3D 
reconstruction model is updated according to their difference. 

The forward projection to a detector pixel i at a projection angle can be used to illustrate 
15 the whole forward projection problem. An x-ray beam containing N i photons is incident from 
the source to the center of the selected detector pixel. This beam penetrates a series of object 
voxels and is sequentially attenuated by them. The total aggregate attenuation is < l,u (n) > i and 

the number of transmitted photons is N | .e~ </ '" ( " )> ' , which is the forward projection to the pixel. 

This operation is repeated for all detector pixels that form the forward projection at this angle. 
20 The forward projections at all angles can be done in the same way except that the "pseudo-beam" 
is rotated. 

The 3D reconstruction model is updated at the backprojection step. Equation 3 describes 
the update of a voxel j at the n-th iteration of reconstruction. The whole image is updated by 
25 doing the same operation on every voxel in it. At a projection angle, the center of the voxel is 

projected from the source to a detector pixel containing the attenuation information of this voxel. 
This operation is repeated at other angles and totally 1 1 detector pixels are found. In equation 3, 
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the values of these 1 1 pixels, both in forward projection and in measured projection are used to 
update the object (reconstruction) voxel (the summation is on the set of these 1 1 pixels). 

3 J Positions of x-ray source, object voxel and detector pixel 
5 The origin of the coordinate system is at the axis of rotation 14 as illustrated in Figure 2. 

The rotation plane of the x-ray 12 source is the YZ-plane (x=0). The detector 20 is parallel to the 
XY-plane at z=21 .7cm. The distance between the source 12 and the axis of rotation 14 is D sa and 
the distance between the detector 20 and the axis of rotation 14 is D da . At projection angle 9, the 
position of the x-ray source 12 is: 

*,(0) = o 

10 y s (e) = D sa sm(9) (4) 

Z,{0) = -D M .GQS(e) 

The reconstructed object 24 is a rectangular volume, represented by a three-dimensional 
array of voxels 26. The breast volume 18 is contained in this rectangular volume 24. In a 
reconstructed image, the value of a voxel is positive if it represents breast tissue; zero if it 
represents the empty space out of the breast. In the coordinate system, the position of a voxel 26 
15 indexed by (m x , m y , m 2 ) is: 

. X obj= X obj+ m x' d x 

yobj= Y o bj +™ y -d y (5) 

Z obj= Z obj+ m z- d z 

where {X 0 bj, Y 0 bj, Z 0 bj) is the position of the center of the rectangular volume 24; d X9 d y and d z are 
the size of the voxel 26 in three dimensions. 

20 The position of a detector pixel 28 indexed by (n x , n y ) is: 

X p= X p+ n x' d 'x 

yp = Y P+ n y-<> ( 6 ) 

J 

where (X P9 Y p ) is the position of the center of the detector 20; d\ and d y are the size of the pixel 
28 in X and Y dimensions. 
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3.2 Forward projection 

The forward projection is implemented by ray tracing from the x-ray source 12 to 
detector pixel 28. At a projection angle, the x-ray beam to a detector pixel 12 is attenuated from 
5 the point where the beam enters the volume 24 to the point where it goes out. The total 

attenuation along the beam < /, >. is calculated by accumulating the attenuation / • u (n) by 

each voxel 26 on the beam line. The number of transmitted x-rays to the pixel 28 is A/ r | .e" </,M(B>> ' . 

The forward projection of the object 18 at this angle is obtained by repeating this operation for 
all detector pixels 28. The forward projections at other angles are calculated in the same way 
10 except the x-ray source 12 is at a different location. 

The first step of forward projection is to determine the orientation of the x-ray beam 30 as 
illustrated in Figure 3. At an angle, the position of the x-ray source (x s ,y Sy z s ) 12 and detector 

pixel (x p , y p , z p ) 28 are determined by equation 4 and 6. The orientation of the beam P(x, y ,z) 30 
15 from source 12 to the detector pixel 28 can be described by two parameters: (1) P, the angle 

made by the beam and the YZ-plane; (2) a, the angle made by the projection of the beam in YZ- 
plane and the Z-axis. These two parameters are determined by: 

a = tan- x ((y p -y s )/(z p -z 5 )) 

(7) 



P = tan" 1 ^ /j(y p -y s y+( Zp - Zs y) 



The path length P X P 4 of the x-ray beam 30 through a voxel 26, as illustrated in Figure 4, 
20 is also the distance between the centers of two successive voxels along the beam. The position 



of the next voxel along the beam can be located by shifting Ax , Ay and Az (P X P 2 , P 2 P 3 and 



P 3 P 4 in Figure 4) along three dimensions from the current voxel 26. 



Ax = P X P 2 = i^P 4 .cos p.cosa 



Ay = P 2 P 3 =P ] P 4 .cosj3.sina (8) 



Az = P 3 P 4 = P l P 4 .sia/3 
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To calculate 
calculated first: 



P^P 4 , its projection in the YZ-plane P X P 3 , illustrated in Figure 5, is 



PP 

13 



= d y / since if a > tan \d y /d 2 ) 9 



PP 



= d 2 Icosa if a < tan /d 2 ) 



(9) 



In a similar way, the path length P l P A can be calculated by: 



P t P 4 = d x I sin J3 if y? > tan"' / P,P 3 
^U^/cos^ if ^^tan _l (rf,/^P 3 



); 
) 



(10) 



10 There are exceptions to the two cases illustrated in Figure 5. In a case shown in Figure 6, 

the path lengths through voxel 3 and 4 cannot be described by equation 10. But the total path 
length of them is equal to the path length in voxel 2. The total attenuation by voxel 3 and 4 is 
equivalent to the attenuation by the shaded area in Figure 6, which has the same path length as 
voxel 2. The equivalent attenuation is estimated by a linear interpolation of attenuations by 

15 voxel 3 and 4. The weighting for the interpolation is proportional to the inverse of the distance 
from the voxel center to the beam line. The ratio of the weighting for voxel 3 to that for voxel 4 
is d4/d 3 , equivalent to r^, where d 3 and d* are the distances from the voxel center to the beam; 
r 3 and r 4 are the distances from the voxel center to the projection of the beam along the Y-axis. 



20 The total attenuation along a beam to a detector pixel i is the summation from the first 

voxel at the point where the beam enters the volume to the voxel at the point where the beam 
goes out of the volume. For a beam with orientation (a, (5), the position of the voxel at entering 
point is: 

25 y Q =y s +(z 0 -z s )tma; (11) 

z 0 =21.7-2); 
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where D is the thickness of the reconstruction volume. The attenuation / • u 0 by the first voxel at 
(x0, yo, z 0 ) is calculated and then the tracing point is shifted forward by (Ax, Ay, Az) to the next 
voxel along the beam, where the attenuation l>u { is calculated and added to / • u 0 . At the n-th 
step, the position being search is: 

x n = x 0 + n • Ax 

5 y m =y*+n-6y (12) 

z n = y 0 +n-Az 

The number of steps of forward projection is V = int(Z) / Az) + 1 . After V steps, the total 

M 

attenuation along the beam to detector pixel i is ^ u n • l n (represented by < /, u {n) > i ). The 
number of transmitted x-ray photons is N ( e ns0 " " . 
3.3 Backprojection 

The value of the object voxel is updated at the backprojection step as illustrated in Figure 
7. At this step, projection pixels containing the attenuation information of the selected object 
voxel are found and used to update the value of this voxel. At a projection angle, the position of 
the detector pixel (x p , y p , z p ) which contains the information of a selected voxel is: 

* P = *s + (*obj -Xs)< z p- Z s) Kz obJ - Z s ) 

y P = y s + (y ob j -y s )( x P - x s ) - ^ ) (i 3 ) 

z p =2\n 

where (x 5 , z 5 ) is the position of the x-ray source at this angle. This operation is repeated to 
find detector pixels related to this voxel at other angles. The value of this voxel is updated by 
equation 3, using these detector pixels. 

4. Image Reconstruction Results 

4. 1 Study on an ACR Phantom/Mastectomy Specimen 

A phantom 38 is composed of a piece of mastectomy specimen 40 and a feature plate 42 
25 from an American College of Radiology (ACR) accredited mammography phantom and placed 
on detector 20 as illustrated in Figure 8A. The feature plate 42, further illustrated in Figure 8B, 
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contained nylon fibers (labled 1 to 6 on the plate), simulated micro-calcifications (labeled 7 to 1 1 
on the plate) and tumor-like masses (labeled 12 to 16 on the plate). The mastectomy specimen 
40 is a surgically removed breast tissue containing lesions. The combination of the feature plate 
42 with the mastectomy specimen 40 makes it very hard to find features of the ACR phantom 42. 
5 The reconstructed feature plate demonstrates how the three-dimensional reconstruction works to 
improve the visibility of features. 

Ten features (fiber 1, 2, 3, 4; micro-calcification cluster 7, 8, 9 and mass 12, 13, 14) can 
be seen very well in a projection of the 4cm thick ACR phantom 42 itself (Rh/Rh, 28kVp and 
10 160mAs) as shown in Figure 9A. With the superimposed mastectomy specimen 40, only one 
feature (micro-calcification cluster 7) is visible in a projection (Rh/Rh, 30k Vp and 140mAs) as 
can be seen in Figure 9B. 

The reconstruction of the feature layer after 10 iterations is shown in Figure 9C. The x- 
15 ray energy and exposure are the same as that used to create the image of Figure 9B. More 

features (micro-calcification cluster 7, 8, 9 and mass 12) can be seen in the reconstruction. Even 
some low contrast features (fiber 1, 2, 3, 4) are recognizable. The number "503 059" on the label 
is clearer. It is clear that the visibility of features are significantly improved. 

20 4.2 3D Reconstruction of a Patient Tissue 

Clinical imaging of volunteers conducted at Massachusetts General Hospital under IRB 
approved protocols have been reconstructed for comparison of conventional film-screen 
mammography and to tomosynthesis mammography. As an example, a mediolateral oblique 
(MLO) mammogram from a volunteer was obtained using film-screen system (Mo/Mo, 25kV 

25 and 330mrad average glandular dose). The x-ray film image is shown in Figure 10. The patient 
was found to have a non-palpable 10 mm invasive ductal cancer with associated in situ tumor 
and this was proved by biopsy. The cancer was difficult to see in the conventional screening 
mammogram and was found primarily because the calcifications associated with it drew the 
attention of the radiologist. 

30 
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A tomosynthesis image dataset was taken with Rh/Rh target/filter at 28k Vp and a total 
dose of 307mrad. Three reconstructed slices from the 3D reconstruction are shown in Figure 11. 
Blood vessels are seen near the breast skin in Figure 1 1 A. A tumor that has intraductal as well as 
invasive ductal cancer elements is just out of the plane of section in Figure 1 IB. The invasive 
5 tumor mass, marked by an arrow, with associated calcifications in the in situ portion is clearly 
seen in Figure 1 1C, as is a benign intramammary lymph node in the upper portion of the image. 

It is apparent from this volunteer's dataset that overlapping structures in the conventional 
two-dimensional projection images (Figure 10) were spacially separated. A reconstructed image 
10 provided at three different depths (Figure 1 1 A illustrating a depth of Z=2mm, Figure 1 IB 

illustrating a depth of Z=22mm, and Figure 1 1C illustrating a depth of Z=32 mm) makes it easier 
to see the tumor and calcifications and their relative geometry. 

A person of ordinary skill in the art will appreciate further^features and advantages of the 
15 invention based on the above-described embodiments. Accordingly, the invention is not to be 
limited by what has been particularly shown and described, except as indicated by the appended 
claims or those ultimately provided in a utility application claiming priority to this provisional 
application. A number of references have been referred to in the specification by last name of 
the first listed author and year of publication; those references are listed by full citation in the 
20 Bibliography below. All publications and references cited herein are expressly incorporated 
herein by reference in their entirety, in particular, each of the references listed in the 
Bibliography below is expressly incorporated for the teachings referred to in the sections of the 
application above for which they are cited. 



25 



-16- 



Attorney Docket No. 102282-15 



Bibliography 

United States Patent No. 5,872,828 to Niklason et al., entitled "Tomosynthesis System for Breast 
Imaging." 

5 

J. A. Browne, and T. J. Holmes, "Developments with Maximum Likelihood X-ray Computed 
Tomography," IEEE Transactions on Medical Imaging, 1 1(1): 40-52 (1992). 

D. P. Chakraborty, M. V. Yester, G.T. Barnes and A. V. Lakshminarayanan, "Self-masking 
10 subtraction tomosynthesis," Radiology, 150:225-229 (1984). 

P. Haaker, E. Klotz, R. Koppe, R. Linde and H. Moller, "A New Digital Tomosynthesis Method 
With Less Artifacts for Angiography," Medical Physics, 12(4): 431-436 (1985). 

15 T. J. Herbert and R. M. Leahy, "A Generalized EM Algorithm for 3-D Bayesian Reconstruction 
from Poisson Data Using Gibbs Priors," IEEE Transactions on Medical Imaging, 8(2): 194-202 
(1989). 

K. Lange and J. A. Fessler, "Globally Convergent Algorithm for Maximum a Posteriori 
20 Transmission Tomography," IEEE Transactions on Image Processing, 4: 1430-1438 (1995). 

E. Levitan and G. T. Herman, "A Maximum A Posteriori Probability expectation Maximization 
Algorithm or Image Reconstruction in Emission Tomography," IEEE Transactions on Medical 
Imaging, MI-6(3): 185-192 (1987). 

25 

S. H. Manglos, G. M. Gagne, F. D. Thomas and R. Narayanaswamy, "Transmission Maximum- 
Likelihood Reconstruction with Ordered Subsets for Cone Beam CT," Physics in Medicine and 
Biology, 40: 1225-1241 (1995). 

30 L. T. Niklason, B. T. Christian, L. E. Niklason, D. B. Kopans, D. E. Castleberry, B. H. Opsahl- 
Ong, C. E. Landberg, P. J. Slanetz, A. A. Giardino, R. M. Moore, D. Albagi, M. C. DeJule, P. A. 
Fitzgerald, D. F. Fobare, B. W. Giambattista, R. F. Kwasnick, J. Liu, S. J. Lubowski, G. E. 
Possin, J. F. Richotte, C-Y Weinad R. F. Wirth, "Digital Tomosynthesis in Breast Imaging," 
Radiology, 205: 399-406 (1997). 

35 

L. T. Niklason, B. T. Christian, L. E. Niklason, D. B. Kopans, P. J. Slanetz, D. E. Castleberry, B. 
H. Opsahl-Ong, C. E. Landberg, B. W. Giambattista, "Digital Breast Tomosynthesis: Potentially 
a New Method for Breast Cancer Screening," Digital Mammography, edited by N. Karssemeijer, 
M. Thijssen, J. Hendriks and L. van Erning, 13: 51-56 (Kluwer Academic Publishers, 1998). 

40 

T. Pan, B. M. W. Tsui and C. L. Byrne, "Choice of Initial Conditions in the ML Reconstruction 
of Fan-Beam Transmission with Truncated Projection Data," IEEE Transactions on Medical 
Imaging, 16(4): 426-438 (1997). 



-17- 



Attorney Docket No. 102282-15 



A. J. Rockmore and A. Macovski, "A Maximum Likelihood Approach to Transmission Image 
Reconstruction From Projections," IEEE Transactions on Nuclear Science, 24: 1929-1935 
(1977). 

L. A. Shepp and Y. Vardi, "Maximum Likelihood Reconstruction for Emission Tomography," 
IEEE Transactions on Medical Imaging, MI-1: 1 13-122 (1982). 

S. Suryanarayanan, A. Karellas, S. Vedantham, S. J. Glick, C. J. D'Orsi, S. P. Baker and R. L. 
Webber, "Comparison of Tomosynthesis Methods Used with Digital Mammography," Academic 
Radiology, 7:1085-1097, (2000). 

Z. Zhou, R. M. Leahy and J. Qi, "Approximate Maximum Likelihood Hyperparameter 
Estimation for Gibbs Priors," IEEE Transactions on Image Processing, 6(6): 844-861 (1997). 



What is claimed is: 



-18- 



